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I.  Introduction 


The  study  of  incompressible  curved  channel  flow  can  illuminate  some  important  issues  deal¬ 
ing  with  laminar-turbulent  transition.  The  streamwise  curvature  of  the  flow  induces  an 
instability  resulting  in  longitudinal  vortices  (hereafter  referred  to  as  Dean  vortices  [1]).  For 
the  case  in  which  the  gap  between  the  walls  is  small  relative  to  the  radius  of  curvature,  Reid 
[2]  described  two  approximate  methods  which  both  shewed  that  the  vortices  first  become 
unstable  when  the  Dean  number  (De  =  2R eyjd*/r*,  where  Re  =  U*h* /u  is  the  Reynolds 
number,  U*  is  the  bulk  velocity,  d *  is  the  channel  width,  h*  —  d*  j 2  is  the  channel  half-width, 
r*  is  the  inner  wall  radius,  and  v  is  the  kinematic  viscosity)  exceeds  a  value  of  about  36. 
These  vortical  structures  have  been  found  experimentally  by  Ligrani  and  Niver  [3].  The 
linear  and  nonlinear  properties  of  these  axisymmetric  vortices  as  well  as  their  transitions 
into  two  distinct  types  of  wavy  vortices  have  been  well  documented  in  the  direct  numeri¬ 
cal  simulations  and  weakly  nonlinear  studies  of  Finlay,  Keller,  and  Ferziger  [4],  For  more 
mildly  curved  channels  than  they  studied,  the  minimum  Reynolds  number  for  instability  of 
the  Dean  vortices  is  comparable  with  that  for  two-dimensional  (2D)  Tollmien-Schlichting 
(TS)  waves.  In  this  regime,  Gibson  and  Cook  [5]  showed  that  oblique  waves  are  never  the 
dominant  linear  disturbance.  Hence  a  study  of  the  nonlinear  interactions  between  2D  TS 
waves  and  the  streamwise  oriented  Dean  vortices  can  be  meaningfully  conducted. 

Daudpota,  Hall  and  Zang  [6]  (hereafter  referred  to  as  DHZ)  developed  a  weakly  nonlinear 
interaction  theory  to  study  the  interaction  of  Dean  vortices  and  TS  waves  in  curved  chan¬ 
nels.  They  employed  a  multiple  scale  version  of  the  Stuart  [7]  and  Watson  [8]  approach  to 
derive  two  coupled  Landau  equations  for  the  perturbation  amplitudes  of  Dean  vortices  and 
TS  waves.  A  comparison  of  their  theory’s  predictions  with  the  results  of  direct  numerical 
simulation  [9]  suggests  that  their  results  are  in  error  with  respect  to  the  influence  of  the 
TS  wave  on  the  Dean  vortex.  The  resolution  of  this  discrepancy  is  the  major  motivation 
for  undertaking  this  work.  In  section  II.  we  present  a  slightly  different  formulation  of  a 
weakly  nonlinear  interaction  theory  and  describe  how  the  resulting  equations  are  solved. 
The  current  approach  is  based  on  the  work  by  Herbert  [10,  11].  We  compare  some  aspects 
of  this  approach  with  that  used  by  DHZ.  In  section  III.  we  report  some  results  of  the  inter¬ 
action  theory  while  in  section  IV.  we  compare  specific  cases  with  results  of  direct  numerical 
simulation.  Finally  in  section  V.  we  draw  conclusions. 


II.  Mathematical  Formulation 


A.  The  Basic  Equations 


The  incompressible  Navier-Stokes  equations  in  cylindrical  coordinates  (r*,0,  z*)  are  written 
as: 
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and  p*  and  i/*  represent  the  constant  density  and  kinematic  viscosity  respectively.  The 
asterisks  indicate  dimensional  quantities.  The  geometry  of  the  problem  of  azimuthal  flow 
between  infinite  concentric  walls  of  outer  radius  r*  and  inner  radius  r*  is  illustrated  in  Figure 
1.  The  wall  boundary  conditions  require  that: 

U/  =  Us *  =  UT  =  0,  at  r*  =  r*,  (7) 

A  solution  to  the  equations  gives  (Ur* ,Ug* ,UZ* ,  P*)  as: 

(Ur'\  (  _0  \ 

tV  -  (R) 

u:  -  _  o  •  (8) 

\  P*  /  \0p*U*2  J  f(r)2 /r dr  ) 

Here  77*  is  the  bulk  velocity  and 

/(r‘)  =  D*  (r*  log r*  +  <7V  +  P*/r*) 

where: 

c*  =  (r*2  log r*  -  r*2  logr^)  /  (r*2  -  r*2)  (9) 

£*  =  -(ryo)2/K2-r:2)\  ogv  (10) 

D *  =  — 2/(r*  +  E* /h*  log 77) .  (11) 


The  quantity  rc  (rg  -j-  T^'jJ 2  is  the  centerline  radius  and  the  ratio  77  —  r* jr*a  describes  the 
channel  curvature. 

We  nondimensionalize  all  spatial  coordinates  with  the  channel  half-width,  velocities  with 
the  bulk  velocity,  and  pressure  with  p*TF2.  The  temporal  scale  is  h*/TF.  The  Reynolds 
number  is  defined  as 

Re  =  U*h* /u*. 

The  nondimensional  equations  are  easily  obtained  from  the  dimensional  ones  by  replacing  all 
starred  quantities  with  their  corresponding  nondimensional  unstarred  ones  and  noting  that 
p  =  1,  h  —  1,  and  v  =  1/Re.  The  expression  1  -  77  appears  often  below  and  is  denoted  by  A. 
When  A  =  Ac  =  2.179  x  10~5  the  minimum  Reynolds  number  for  instability  of  TS  waves  is 
the  same  as  that  for  Dean  vortices  (see  Gibson  and  Cook[5]). 

We  follow  closely  the  perturbation  method  introduced  by  Herbert  [10,  11],  extending  it 
to  the  case  of  interacting  disturbances.  The  steady  solution  is  perturbed  such  that 

/  Ur  \  /  UT  \ 

Ue  /(r)  -f  ug 

U,  uz  (12) 
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where  ur,  ue,  uz,  and  p  are  perturbation  quantities.  Substituting  the  perturbations  into  the 
nondimensionalized  Navier-Stokes  equations,  subtracting  the  steady  flow  component,  and 
rearranging  gives: 
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where: 
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rdr[r  dr)  +  r2d62  +  dz2  (18) 

The  entire  flow  is  assumed  to  be  periodic  in  the  azimuthal  (streamwise)  and  axial  (spanwise) 
directions;  the  solution  evolves  in  time.  In  addition  the  total  pressure  gradient  is  kept 
constant. 
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B.  Fourier  Expansions 


The  linearized  disturbance  equations  are  obtained  by  setting  the  right  hand  sides  of  Equa¬ 
tions  13-16  equal  to  zero.  A  solution  to  the  equations  for  a  2D  TS  wave  can  be  written 


«r1,0(r»*) 

U81,0(r,t) 

0 


exp  i(aQ  —  7 (t)). 


Here  a  is  real  and  represents  the  azimuthal  wavenumber  while  j(t)  =  got  is  a  real  function 
defining  the  phase.  Any  growth  or  decay  of  the  TS  wave  is  expressed  in  terms  of  an  amplitude 
parameter  such  that: 


uV  -°(rti) 

0 

P1,0(r)  t) 


=  A(t ) 


*  1,0, 0,0/ \ 

Ug  (.r  J 


p1  ,°,o,0(r) 


where 


dA 

li  ~  a°'°-4 


and  a0,o  is  the  linear  growth  rate.  The  apparently  superfluous  superscripts  and  subscripts 
are  included  here  for  consistency  with  the  notation  used  later.  A  convenient  normalization 
for  this  problem  is 

|wr1,0(rc,  i)i  =  A(t)  (22) 

where  rc  is  the  channel  centerline.  We  take 

ur1M0(re)  =  1,  (23) 

thus  fixing  the  amplitude  and  phase  of  the  TS  wave. 

In  a  similar  manner,  one  can  write  the  solution  for  the  primary  Dean  vortex  mode  as: 

/ur\  /  dr0,1(r,  t)  \ 

u*  =  !z expf(/?z-C(0)»  (24) 

ut  uz  ( r,t ) 

\  p  J  \  t )  / 

where  (3  is  the  real  axial  wavenumber  and  £( t )  =  h0t  =  0  indicates  that  the  linear  Dean 
disturbance  is  steady.  We  also  write 
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where 


(26) 


dB 

Ti  =  b°-°B- 


The  normalization  is  chosen  such  that: 


\Ur°''('c,t)\  =  B(t), 


(27) 


or  in  this  case, 

*w*(re)  =  1.  (28) 

It  is  conceivable  that  a  particularly  inauspicious  choice  of  the  parameters  Re  ,  A,  and 
wavenumbers  might  lead  to  uTl'0'0fi(rc)  =  0  or  ur0,0,1,0(rc)  =  0  and  hence  the  normaliza¬ 
tion  used  in  either  Equation  22  or  27  could  fail.  Such  an  event  would  lead  to  overflow  errors 
in  the  numerical  procedure  and  an  alternative  normalization  could  be  used.  Herbert  [11] 
discusses  integral  normalizations  which  are  more  robust  but  somewhat  more  complicated  to 
implement.  In  practice  we  have  never  experienced  any  problems  with  the  normalizations 
in  Equations  22  and  27.  We  attribute  this  to  the  fact  that  the  amplitudes  of  the  radial 
velocities  of  the  dominant  TS  waves  and  Dean  vortices  do  not  have  local  minima  near  the 
channel  center. 


When  the  nonlinear  terms  (the  right  hand  sides  of  Equations  13-16)  are  included,  the 
disturbances  can  interact  with  themselves,  the  mean  flow,  each  other  and  all  their  relevant 
complex  conjugates.  This  results  in  the  generation  of  harmonics,  mean  flow  distortions, 
and  various  corrections  to  the  fundamental  disturbances.  Hence  it  is  natural  to  expand  the 
perturbation  variables  in  a  double  Fourier  series: 


(uT\ 

Ug 

Uz 

\  P  ! 


E  E 

n=  — oo  I  =  — oo 


(Urn’\r,t)\ 
Ugn'l(r,  t ) 
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exp  in(aff  —  j(t))  exp  il(/3z  —  ((t)). 


(29) 


C.  Amplitude  Expansion 


The  system  of  nonlinear  partial  differential  equations  obtained  by  substituting  Equation  29 
into  Equations  13-  16  are  coupled  and  difficult  to  solve  efficiently.  We  seek  a  solution  by 
expanding  in  the  amplitude  parameters  A(t )  and  B(t)  about  the  linear  solutions  given  by 
Equations  19  and  24.  The  solutions  will  reflect  the  deviations  from  linear  behavior  for  finite, 
but  sufficiently  small  amplitudes  A  and  B.  In  this  context,  all  harmonics  of  the  primary 
disturbances  are  considered  to  be  forced.  Since  the  linear  TS  and  Dean  solutions  are  0(A) 
and  0(B)  respectively,  their  first  harmonics  and  cross  terms  will  be  0(A2),  0(B2),  and 
O(AB).  Higher  order  harmonics  and  cross  terms  will  be  0(AnBl),  hence  it  is  reasonable  to 


let: 


/  uTn’l(r,t)  \ 
uen'l(r,  t ) 
uzn'l(r,  t ) 
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A\n\BV‘ 
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(30) 
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where  all  double  tilde  terms  are  0(1)  except  for  those  with  n  =  l  =  0.  Heuristically,  one 
can  see  that  the  exception  with  the  n  =  l  =  0  term  comes  about  from  the  fact  that  the 
lowest  order  mean  flow  distortions  are  generated  from  the  product  of  either  TS  or  Dean 
fundamental  disturbances  with  their  respective  complex  conjugates.  Since  the  fundamental 
disturbances  are  0(A)  and  0(B) ,  the  product  terms  which  generate  n  —  l  =  0  are  0(A2) 
and  0(B2)  respectively.  Substitution  of  the  representation  in  Equation  30  into  Equation  29 
gives: 


-EE 

n=-oo  l=- oo 


urn'l(r,t) 
s  n.!/„ 


Uzn,if’tl  A}n\t)B',l\t)  exp  in(a&  —  j(t))  expil((3z  —  £(£)).  (31) 

u2  (r,t) 


~n’h  \ 
p  (M) 


Since  we  seek  solutions  which  are  small  deviations  from  the  linear  behavior,  we  now  expand 
the  double  tilde  representations  in  sums  which  are  the  products  of  ascending  powers  of  the 
amplitude  functions  with  coefficients  which  are  strictly  functions  of  the  radial  coordinate. 
In  the  limit  as  A  — *  0  and  B  -*  0,  the  solutions  tend  towards  the  linear  results.  Because  of 
the  invariance  of  the  original  equations  and  boundary  conditions  with  respect  to  arbitrary 
translation  in  the  streamwise  and  spanwise  directions  and  the  assumption  of  periodicity  of 
the  solutions  in  these  directions,  only  even  powers  of  A  and  B  are  needed  [12].  Hence, 

'  irn\r,  t)\  /  urn'm'l'k(r)  \ 

“Vi  '  =EE  “V.J  (32) 

( r,t )  m=0k=0  u* 

p‘(r,t)j  P"UW  ' 

By  substitution  into  Equation  31  one  obtains  the  full  representation  of  the  perturbation 
variables: 


I  ur  '  '  [r ) 

oo  oo  oo  oo  o  Ti.m.l.kf  \ 

_  V'  V  Ue 

ZZ  Z-,  Z-,  Z^  n,m,l,k/  \ 

n—  oo  m~0 /=  — oo  fc=0  z  '  ' 


A2m+]n](t)B2k+W(t) 


exp in(a6  —  7 (t))  expi l((3z  —  £(t)).  (33) 

In  order  that  the  wall  boundary  conditions  be  satisfied  at  all  orders  of  approximation,  we 


require: 


\  (0 

uen-m-l-k(r)  =  0 

uzn'm’l'k(r)  J  VO 


at  r  =  r„  r0 


for  all  combinations  of  n,  m,  /,  and  k. 


The  time  derivatives  A/ A,  B  / B,  7,  and  (  are  also  expanded  in  terms  of  A2  and  B2  such 


YJY,arn,kA2mB2k 


771=  0  A:— 0 
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(35) 


D  OO  OO 

r  =  J2Eb^A2mB,k  M 

0  m=0  k= 0 

7  =  f^f^gmikA2mB2k  (37) 

m=0  k~0 

OO  OO 

C  =  EEM2m52k  (38) 

771  =  0  fc=0 

Equations  35  and  36  are  the  coupled  Landau  equations  which  describe  the  growth  of  the 
disturbances.  The  purpose  of  the  remainder  of  this  section  is  to  determine  the  Landau 
coefficients  am  fc  and 


D.  Solution  Method 


The  solution  expansions  in  Equations  31-38  are  substituted  into  the  perturbation  equations 
13-16.  All  terms  with  common  exponential  factors  are  grouped  together.  The  simplest 
nonlinear  theory  is  obtained  by  considering  only  those  terms  with  Fourier  exponents  limited 
by: 

—2<n,l<2. 

With  this  restriction,  it  is  only  appropriate  to  consider  terms  in  the  amplitude  expansions 
(Equation  32)  with  m  =  0,1  and  k  =  0,1.  Larger  values  of  m  or  k  give  higher  order 
contributions  which  should  not  be  considered  at  this  degree  of  truncation.  In  each  group 
with  common  exponential  factors,  all  terms  with  common  powers  of  A  and  B  are  collected. 
This  tedious  process  of  expanding  the  solutions,  substituting  into  the  equations  and  collecting 
terms  with  common  factors  of  exponentials  and  amplitude  parameters  is  done  automatically 
using  the  symbolic  manipulator  Macsyma.  The  case  n  =  m  =  1  =  fc  =  0  just  represents  the 
basic  flow  which  has  already  been  subtracted  out.  A  sequence  of  12  systems  of  equations 
remains.  Each  system  contains  four  equations  and  four  unknown  expansion  coefficients 
which  are  functions  of  r.  A  description  of  each  of  the  systems  is  summarized  in  Table  1. 
The  equation  system  number  indicates  the  order  in  which  we  chose  to  solve  the  equations. 
To  ensure  that  all  information  needed  to  solve  any  system  has  already  been  computed,  all 
systems  of  lower  order  in  the  amplitude  parameters  should  be  solved  before  those  of  higher 
order. 

The  first  two  equation  systems  are  eigenproblems  of  the  form: 

Ln,/Qn,o,/,o  =  +  (^(oo.o  —  f<7o,o)  +  ^  (^o,o  ~  ih o,o))  Mq^o.j.o  (39) 


where: 


f  Ur  (7) 

Ug  (r) 

w*  (r) 

pnmJ.t(r) 


(40) 
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Eq.  Sys.  No. 

Order 

n 

m 

1 

X 

Description 

1 

A 

1 

0 

0 

0 

Linear  TS  wave 

2 

B 

0 

0 

1 

0 

Linear  Dean  vortex 

3 

A2 

0 

1 

0 

0 

Mean  flow  distortion  from  TS 

4 

B2 

0 

0 

0 

1 

Mean  flow  distortion  from  Dean 

5 

A2 

2 

0 

0 

0 

TS  harmonic 

6 

B2 

0 

0 

2 

0 

Dean  harmonic 

7 

AB 

1 

0 

1 

0 

TS  -  Dean  cross  harmonic 

8 

AB 

1 

0 

-1 

0 

TS  -  Dean  cross  harmonic 

9 

A3 

1 

1 

0 

0 

Self-correction  to  TS 

10 

B3 

0 

0 

1 

1 

Self-correction  to  Dean 

11 

AB2 

1 

0 

0 

1 

Dean  correction  to  TS 

12 

A2  B 

0 

1 

1 

0 

TS  correction  to  Dean 

Table  1:  Summary  of  equation  systems 


and 


f-l 

0 

0 

\  o 


0  0  0\ 

-10  0 
0-10 
0  0  0/ 


(41) 


For  the  eigenproblems,  the  indices  n  and  l  may  only  have  the  values  0  and  1  with  the 
constraint  n  +  l  =  1.  The  elements  of  the  operator  LnI  are  given  in  Appendix  A.  The 
solution  of  the  eigenproblem  determines  the  eigenvalues  (we  consider  here  only  the  least 
stable  ones),  a0,o  +  Afo.o  and  b0to  +  ih0lo-  The  eigenvectors  qi ,o,o,o(’')  and  qo,o,i,o(n)  are  known 
to  within  a  constant  factor.  The  normalizations  in  Equations  23  and  28  are  used  for  the  two 
eigenvectors.  The  real  parts  of  the  eigenvalues  are  the  zeroth  order  Landau  coefficients. 


Equation  systems  3  and  4  represent  the  mean  flow  distortion  caused  by  the  TS  and 
Dean  disturbances  respectively.  For  this  case,  n  =  l  =  0  and  the  operator  Lnii  can  be  greatly 
simplified.  Here  we  consider  only  the  case  where  either  m  or  k  equals  1,  while  the  other  equals 
0.  Earlier  we  specified  that  the  mean  pressure  gradient  is  constant,  p0,1,0,0  =  p0’0,0,1  =  0. 
Analysis  indicates  that  ur0,1,0’ '  =  iq.0,0,0,1  =  iq0,1,0,0  =  iq0,0,0,1  =  0.  The  remaining  equation 
can  be  written  as: 


1  d  *  0,m,0,k  1  <P  *  O.m.O.fc 

“  Re*7  5 

-  Ug0’0,1’0!/*0’0'-1’0) 

,7  m,0,fc,0  -  -m,0,-fc,0\ 


(2mao,o  +  2A:i0,o  +  r'2j^e)^0’m’°’fc 

ikp{ue0-°'-l’0uz0M 

-  (ug~m'0'-k’0uTTU'0'k'°  +  uem-0'k'°ur-m'0’-k'Q)/r 
-  (£ (ue-mfl'-k’0)urm'°'k'0  +  ±(uem'0'k'Q)ur-m'0'-k'0).  (42) 

An  interesting  restriction  to  the  theory  results  from  Equation  42.  Consider  the  associated 
homogeneous  problem: 

1  d<f>  1  d?cf) 


(2ma0,o  +  'Zkboo  + 


r2Re 


)<f> 


rRe  dr  Re  dr 2 


=  0 
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with  boundary  conditions: 


(f>  =  0  at  r  =  rt,  ra. 


Apply  the  transformation: 


r2  =  (2 /  (-(2mao,0  +  2&&0,o)Re) 

to  the  homogeneous  problem  and  divide  by  (2mao,o  +  2kb00)  to  obtain: 


k  (dt  d£ 2 


=  0 


with 

where 


4>  =  0  at  £  =  &, 


(43) 


(l  =  r?/  (— (2ma0,o  +  2kb0fi)Re ) 

C  =  rH  (-(2ma0l0  +  2fc&0,0)Re) 


Equation  43  is  the  Bessel  equation  of  order  1.  If  the  points  <f,  and  correspond  to  zeros  of 
the  appropriate  Bessel  function,  then  this  free  eigenmode  can  contribute  an  arbitrarily  large 
component  to  the  solution  of  the  forced  problem.  Moreover,  a  numerically  ill-conditioned 
system  arises  whenever  &  and  are  both  in  the  neighborhood  of  zeros  of  the  Bessel  function. 
A  sure  way  to  prevent  this  problem  is  to  require  both  o0, 0  and  60,0  to  be  non-negative.  This 
essentially  restricts  the  application  of  the  theory  to  those  regimes  in  Reynolds  number  - 
wavenumber  space  where  both  the  TS  and  Dean  disturbances  are  unstable.  In  this  regime, 
solving  for  ue°,1,0’°  and  ue0’0'0,1  is  straightforward. 

The  equation  systems  5-12  are  all  of  the  form: 

(Ln  /  —  (2mao,o  +  Mao.o  —  i™<7o,o  +  2&£>o,o  +  Kl^o.o  —  ^0,0)  f,)t  =  R,,  (44) 

where  R,  is  the  right  hand  side  vector  of  the  5  equation  system.  The  right  hand  side  vectors 
consist  of  nonlinear  combinations  of  the  solutions  obtained  at  lower  order.  The  elements  of 
the  right  hand  side  vectors  are  given  explicitly  in  the  Appendix  B.  Some  important  consid¬ 
erations  in  solving  the  equations  are  given  below.  (Many  of  these  issues  are  addressed  in 
detail  by  Herbert  [11]  for  the  case  of  a  single  2D  TS  wave  in  a  straight  channel.) 

With  the  restriction  that  do.o  and  b 0,0  be  non-negative,  the  homogeneous  problems  asso¬ 
ciated  with  Equation  44  are  unlikely  to  have  free  modes.  The  occurrence  of  any  free  modes 
for  Dean  harmonics  would  require  that  {2k  +  l)b0< 0  be  an  eigenvalue  for  the  homogeneous 
problem  with  wavenumber  1/3.  (Note  that  h00  =  0.)  Such  an  unlikely  coincidence  would 
lead  to  obviously  nonsensical  results  of  the  numerical  solution.  TS  harmonics  and  crossed 
TS-Dean  harmonics  would  be  even  less  likely  to  have  free  modes  because  both  real  and 
imaginary  parts  of  the  eigenvalue  would  need  to  match. 


9 


Equation  systems  5-8  for  the  various  harmonics  require  numerically  solving  a  sequence 
of  linear  ordinary  differential  equation  systems.  In  the  case  of  Dean  disturbances,  it  is 
possible  that  the  harmonic  of  an  unstable  mode  is  also  in  the  unstable  regime.  When 
this  happens,  it  violates  the  assumption  that  all  of  the  higher  harmonics  are  forced  by  the 
primary  disturbances.  While  numerical  results  can  be  obtained  for  such  cases,  the  results 
are  meaningless  since  the  harmonic  can  grow  on  its  own  and  is  not  simply  forced  by  the 
fundamental. 

Equation  systems  9-12  define  the  corrections  to  the  fundamental  TS  and  Dean  distur¬ 
bances  and  allow  us  to  calculate  the  Landau  coefficients  ai|0,  &o,i>  ao,i,  and  respectively. 
It  is  in  the  solution  of  these  terms  that  the  method  devexoped  by  Herbert  [10,  11]  differs 
significantly  from  the  standard  Stuart  [7]  and  Watson  [8]  approach. 

The  right  hand  side  vectors  R,  for  s  =  9-12  can  be  rewritten  as. 

R,  =  R'  +  /c.Mq„tcuo  (45) 


where: 


«9  =  <*i,o  —  *01,0 

^10  =  &0,1  *fco,i 

K-ll  —  <*0,1  —  *00,1 

/Cl 2  =  bi,0  -  *Al,0- 


The  correspondence  between  the  index  s  and  the  indices  n,  m,  l,  and  k  is  given  in  Table  1. 
Note  that  these  equations  have  either  n  =  1,  /  =  0,  or  n  =  0,  l  —  1.  No  other  combinations 
of  n  and  l  appear.  The  solution  vectors  qn,w,*  for  these  equation  systems  are  the  sums  of 
particular  solution  vectors,  hence  it  is  useful  to  write: 


£[n,m,l,k{r)  ~  KsXn,m,l,k{r )  Xn,m,l,k{r ) 


(46) 


where: 

Xn,m,i,Jt(r»)  =  Xn,m,l,k{ro )  ~  Xn,m,J,fe(7’*)  —  Xn,m,l,k{ro )  —  0- 

Substitution  of  Equations  45  and  46  into  44  leads  to: 


(47) 


L „i  -  (2ma0,o  +  n(a0i0  —  *0o,o)  +  2fc60,o  +  K^o.o  ~  *^o,o))  M)  —  /CjNlq  n,0,/,0  (48) 

(Ln>j  -  (2ma0,o  +  **(<*o,o  —  *0o,o)  +  2fc50,o  +  ~  *^o,o))  M)  Xn,m,/,fe  =  (4^) 


Recognizing  that 

(Lni  —  (n(a0io  —  *0o, a)  +  l{bo,o  ~  *^o,o))  Nl)  qn.o./.o  =  0 
when  n  =  1,  l  =  0  orn  =  0,  l  =  1,  one  finds  by  inspection  of  Equation  48: 

—  (2ma0,o  +  2A:60,o)  Xn,m,i,fc  =  9n,o,i,o 


(50) 


(51) 
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or 

—  — Qn,o,/,o/  (2ma0io  +  2fci>o,o)  •  (52) 

It  is  straightforward  to  solve  Equation  49  to  obtain  Xn,m  i,k •  Substituting  Equation  52  into 
Equation  46  gives: 

q n,m1/1fc(f)  =  Xn,m,/,fe(r)  ~  *»/  (2mao,o  +  2fc60,o)  qn,o,i,o(r).  (53) 


In  order  to  determine  the  Landau  coefficients  k ,  we  use  the  normalizations  in  Equations 
22  and  27.  Together  with  the  expansions  in  Equations  30  -  33  they  require  that: 


i«V'°(roO!  = 

|  £  jr  A2m+1  =  A{t) 

(54) 

O 

II 

Ji 

o 

II 

6 

I^W)!  = 

|  f)  f^A2mB2k+1ur°^k(rc)\  =  B(t). 

m=0  k= 0 

(55) 

We  already  used  Equations  23  and  28  to  set 

tir  1'0'°’°(rc)  =  1 


and 

Ur0fiX0(rc)  =  1 

hence  for  the  normalization  to  be  correct  at  all  orders  of  approximation, 


tfri,m,°,fc(rc)  =  0 

for  all  m,k  >  0 

(56) 

and 

dr°'m'1’fc(rc)  =  0 

for  all  m,k  >  0. 

(57) 

The  conditions  from  Equations  56  and  57 

can  now  be  used  to  determine  k}. 

Let  the 

components  of  vector  xh,m,i,k  written  as: 


X 


i 

n,m,ltk 


(Wr  \ 

zve 

zuz 

V  ?  / 


(58) 


The  conditions  in  Equations  56  and  57  are  applied  in  Equation  53  to  require: 


t*7,.(rc)  -  Kf/(2ma0,o  +  2kb0,o)  urn'0ll'°(rc)  =  0. 


(59) 


Since  urn'°'l'°(rc)  =  1, 

k,  =  zar(rc)  (2ma0,0  +  2fc60,o)  •  (60) 

Hence  the  Landau  coefficients  k,  are  determined.  The  correction  to  the  fundamentals  can 
be  found  by  substitution  back  into  Equation  46. 
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E.  Numerical  Considerations 


All  of  the  equation  systems  are  solved  using  a  Chebyshev  collocation  procedure  [13].  The 
left-hand  side  operators  for  all  of  the  systems  (except  the  mean  flow  corrections)  are  similar 
and  are  generated  with  a  single  subroutine.  The  operators  for  the  mean  flow  corrections  need 
to  ensure  that  only  streamwise  velocity  perturbations  are  nonzero.  This  was  most  easily  done 
using  a  separate  subroutine.  The  right-hand  side  vectors  were  sufficiently  different  from  each 
other  that  it  was  more  convenient  to  use  Macsyma  to  write  individual  subroutines  for  each  of 
the  equation  systems.  The  two  eigenproblems  are  solved  first  with  a  global  solver  to  obtain 
the  eigenvalue  spectrum.  The  eigenvectors  associated  with  the  most  unstable  eigenvalue  are 
obtained  with  a  local  iterative  procedure.  The  mean  flow  needs  to  be  computed  using  128 
bit  accuracy  for  the  very  weak  curvatures  that  we  used.  We  found  that  65  collocation  points 
across  the  channel  (including  the  endpoints)  produced  Landau  coefficients  which  agreed  to 
at  least  3  digits  with  those  obtained  with  97  collocation  points.  All  results  presented  here 
were  obtained  using  65  points. 


F.  Comparison  of  Approaches 

Here  we  compare  some  aspects  of  the  current  approach  with  that  of  DHZ.  Since  they  do  not 
use  the  normalizations  in  Equations  22  and  27,  they  do  not  have  the  conditions  of  Equations 
56  and  57  and  hence  need  some  alternative  constraint  to  obtain  Landau  coefficients.  They 
use  orthogonality  conditions  to  obtain  the  Landau  coefficients  when  both  Dean  and  TS 
disturbances  are  on  their  respective  neutral  stability  curves.  Near  the  neutral  stability 
curves,  the  linear  growth  rates  vary  linearly  with  Reynolds  number  deviation  off  of  the 
curves.  Deviations  from  the  critical  curvature  parameter,  Ac,  are  also  assumed  to  affect  the 
linear  growth  rates  linearly.  They  find  that  the  effect  on  the  coefficients  of  the  nonlinear 
terms  in  the  Landau  equations  are  independent  of  deviations  away  from  the  neutral  Reynolds 
numbers  and  critical  curvature  to  the  order  considered. 

While  the  two  approaches  involve  the  same  type  of  expansion,  the  different  ways  of  solving 
for  the  Landau  coefficients  makes  each  appealing  for  obtaining  certain  types  of  information. 
Below  we  discuss  some  of  the  restrictions  and  advantages  of  each  method. 

Both  theories  assume  that  the  fundamental  disturbances  are  the  major  structures  in  the 
flow  and  that  none  of  their  harmonics  or  product  wavevectors  are  linearly  unstable.  These 
secondary  waves  must  be  stable  even  in  the  presence  of  finite  amplitude  primary  waves. 

Both  expansions  are  only  valid  for  small  amplitude  disturbances.  Exactly  how  small  is 
not  known  a  priori.  In  order  for  the  expansions  to  be  valid  for  large  times,  it  is  necessary 
for  the  linear  growth  rates  to  be  “small”.  Small  linear  growth  rates  exist  near  the  neutral 
stability  curves.  In  this  sense,  the  restriction  of  being  near  the  neutral  stability  curves  with 
the  DIIZ  approach  is  similar  to  the  requirement  in  the  present  approach  that  the  amplitudes 
be  small.  The  present  approach  might  be  somewhat  more  accurate  if  one  were  not  interested 
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in  long-time  valid  solutions;  e.g.  if  one  wanted  only  to  know  the  first  effect  of  nonlinearity 
on  the  strength  of  an  instability.  In  such  a  case,  the  present  amplitude  expansion  could  be 
valid  in  cases  with  large  linear  growth  rates  as  long  as  the  amplitudes  were  still  sufficiently 
small. 

The  theory  developed  by  DHZ  describes  the  weakly  nonlinear  behavior  of  the  flow  for  a 
range  of  Reynolds  number  and  curvature  ratio  perturbations  away  from  the  neutrally  stable 
conditions.  In  particular,  the  perturbations  could  be  taken  into  the  linearly  stable  regime; 
there  the  Herbert  approach  may  fail  due  to  a  resonance.  The  theory  of  DHZ  does  not  have 
this  problem  since  all  of  their  constants  are  computed  right  on  the  neutral  stability  curves. 
On  the  other  hand,  the  present  theory  may  be  applied  to  any  combination  of  wavenumbers 
and  Reynolds  numbers,  while  the  DHZ  theory  in  its  current  form  works  only  when  the 
azimuthal  and  axial  wavenumbers  correspond  to  neutrally  stable  disturbances  at  the  same 
Reynolds  number. 

The  derivation  and  solution  of  the  equations  of  the  present  theory  are  somewhat  simpler 
than  those  of  DHZ  since  the  present  theory  does  not  require  the  use  of  adjoint  eigenfunctions 
and  the  orthogonality  condition.  Additionally,  because  the  amplitude  function  is  uniquely 
defined  in  terms  of  the  fundamental,  the  change  in  shape  of  the  fundamental  disturbance  can 
be  computed.  This  allows  the  present  theory  to  be  extended  to  higher  orders  unambiguously. 
Herbert  [10]  presents  seven  Landau  constants  for  plane  Poiseuille  flow  using  this  approach. 
With  modified  versions  of  the  same  procedure  used  here,  Singer,  Meyer,  and  Kleiser  [14]  and 
Ng,  Singer,  Henningson,  and  Alfredsson  [15]  compute  multiple  Landau  coefficients  for  the 
development  of  vortices  in  a  three-dimensional  boundary  layer  and  a  rotating  channel  flow, 
respectively.  Unfortunately,  the  derivation  of  the  equations  for  the  higher  Landau  coefficients 
for  the  present  interacting  wave  case  is  too  complicated  to  be  handled  within  the  framework 
of  our  current  Macsyma  program  and  is  not  done  here. 

III.  Results 

Most  of  our  results  are  at  the  critical  curvature  parameter  Ac  since  this  is  where  we  expect 
the  effects  of  the  interaction  to  be  greatest.  Near  the  critical  curvature  parameter,  A  = 
Ac  =  2.179  x  10-5,  we  found  that  small  changes  in  the  curvature  parameter  do  not  change 
the  qualitative  properties  of  the  results.  To  orient  the  reader,  Figures  2  and  3  are  neutral 
stability  curves  for  the  TS  and  Dean  waves  respectively  at  the  critical  curvature.  Regions 
where  a0  0  and  bQ  0  are  positive  represent  the  linearly  unstable  areas.  All  of  our  results  will 
be  for  points  that  are  in  the  unstable  regime  of  both  the  TS  and  Dean  disturbances. 

Figures  4  and  5  illustrate  the  zero  contours  of  the  Landau  coefficients  ali0  and  60,i  in 
regions  where  the  TS  and  Dean  disturbances  are  linearly  unstable,  respectively.  These 
coefficients  represent  the  self-interaction  effects  of  the  TS  and  Dean  disturbances.  Positive 
values  indicate  that  increased  amplitude  exasperates  the  instability.  Negative  values  are 
necessary  for  an  individual  disturbance  to  reach  a  saturation  state  in  the  linearly  unstable 


region.  The  TS  waves  have  negative  self-interaction  coefficients  near  the  lower  branch  of 
the  TS  neutral  stability  '-urve,  but  positive  values  elsewhere.  At  low  Reynolds  numbers,  the 
Dean  disturbances  have  negative  self-interaction  coefficients  for  all  wavenumbers  (3\  at  higher 
Reynolds  numbers  the  self-interaction  coefficient  is  positive  in  the  middle  of  the  unstable 
region. 

Figures  6  and  7  indicate  the  effect  that  interaction  has  on  the  disturbances.  In  both  of 
these  figures,  at  a  given  Reynolds  number,  the  right  and  left  sides  of  the  rectangle  bound 
the  streamwise  wavenumber  range  in  which  individual  2D  TS  waves  are  linearly  unstable. 
Similarly,  for  each  Reynolds  number,  the  top  and  bottom  sides  of  the  rectangle  bound  the 
spanwise  wavenumber  range  in  which  Dean  disturbances  are  linearly  unstable  .  In  Figure  6, 
positive  values  of  a0|1  indicate  that  the  Dean  disturbance  has  a  destabilizing  effect  on  the 
TS  wave.  This  occurs  primarily  in  the  upper  right  hand  corner  of  the  unstable  rectangle, 
though  at  high  Reynolds  numbers,  there  is  also  a  small  region  of  positive  aop  in  the  lower  left 
corner.  Where  a0,i  is  positive  in  the  upper  right  corner,  the  TS  self-interaction  coefficient 
a10  is  also  positive,  so  the  destabilizing  effect  of  a0,i  here  exists  when  the  TS  wave  by 
itself  is  already  quite  unstable.  The  Dean  disturbance  generally  has  a  stabilizing  effect  on 
the  TS  waves  at  low  Reynolds  numbers,  but  as  the  Reynolds  number  increases,  the  Dean 
vortex  has  a  destabilizing  influence  over  a  greater  portion  of  the  domain.  The  corresponding 
results  for  £>10,  the  effect  of  the  TS  wave  on  the  Dean  mode,  are  illustrated  in  Figure  7. 
Low  wavenumber  Dean  disturbances  tend  to  be  destabilized  by  high  wavenumber  TS  waves; 
however  the  fact  that  the  harmonic  of  the  Dean  is  also  unstable  in  these  cases  invalidates 
the  application  of  the  theory  here.  In  no  case  does  the  theory  indicate  that  TS  and  Dean 
disturbances  simultaneously  destabilize  each  other.  Since  the  coefficient,  £>li0  is  positive  over 
progressively  smaller  portions  of  the  domain  with  increasing  Reynolds  number,  it  is  unlikely 
that  simultaneous  destabilization  occurs,  even  at  very  high  Reynolds  numbers.  In  most  of  the 
cases  studied,  the  weakly  nonlinear  interaction  tends  to  stabilize  both  modes.  In  a  physical 
experiment,  multiple  three-dimensional  modes  can  interact,  perhaps  forming  a  secondary 
instability  which  will  eventually  grow  and  lead  to  transition.  A  multi- wave  interaction  and 
subsequent  instability  is  beyond  the  scope  of  the  present  theory. 

For  several  representative  sets  of  parameters,  we  examine  our  results  in  the  context  of 
nonlinear  autonomous  systems.  Recall  that  A  measures  the  amplitude  of  the  TS  wave  and 
D  corresponds  to  the  Dean  disturbance.  Equations  35  and  36  admit  four  possible  steady 
state  solutions:  the  trivial  solution,  finite  A  with  B  —  0,  finite  B  with  A  =  0,  and  a 
combined  state  with  finite  values  of  both  A  and  B.  Each  of  these  states  is  a  critical  point 
and  may  be  classified  according  to  the  local  behavior  of  the  solutions  in  the  vicinity  of  the 
point.  In  the  phase  plane  (A-B),  stable  nodes  are  identified  as  those  for  which  solution 
trajectories  near  the  point  converge  to  the  node,  unstable  nodes  are  critical  points  from 
which  solution  trajectories  diverge,  and  saddle  points  are  those  points  for  which  a  finite 
number  of  trajectories  converge  to  the  point  while  all  others  diverge. 

Several  trajectories  are  illustrated  in  figure  8  for  the  case  where  Re  =  6291.67,  A  = 
2.189  x  10~5,  a  —  74257,  and  f3  =  4.51.  These  parameters  w'ere  used  extensively  by  Singer 
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and  Zang  [9]  in  their  direct  numerical  simulations.  In  the  figure,  the  small  squares  indicate 
the  various  initial  conditions.  An  explicit  fourth  order  Runge-Kutta  differencing  scheme  was 
used  to  trace  the  subsequent  trajectories.  The  filled  circles  show  the  two  stable  equilibrium 
points  at  /I  =  3.8  X  1 0 — 3 ,  B  ~  0,  and  A  =  0,  B  =  7.0  x  10~5  respectively.  Trajectories 
initiated  near  either  of  these  points  tend  to  converge  to  the  respective  nodes.  The  open 
circles  represent  unstable  nodes  and  saddle  points.  Two  initial  conditions  started  quite  close 
to  the  saddle  point  at  A  =  2.4  x  10-4,  B  =  3.5  x  10-5  go  to  different  stable  equilibrium 
states.  T)  e  flow  is  especially  sensitive  to  small  changes  near  the  saddle  point.  The  unstable 
node  at  tue  origin  is  the  endpoint  of  a  semi-infinite  curve  which  goes  through  the  saddle 
point  and  separates  the  region  of  attraction  of  the  two  stable  nodes. 

Figure  9  illustrates  a  typical  trajectory  for  the  case  with  Re  =  4 U 0 0 ,  A  =  2.179  x  10“5. 
a  =  96300,  and  0  =  1.55.  In  this  case  there  are  three  critical  points.  The  origin  and  /I  =  0, 
B  —  4.69  x  10-5  are  both  unstable  nodes.  Trajectories  that  start  in  the  vicinity  of  these 
points  move  away  from  them.  The  point  A  =  3.8  X  10  4,  B  —  6.1  X  10  5  is  a  stable  spiral 
node.  Trajectories  slowly  spiral  towards  this  point.  There  is  no  conventional  equilibrium 
point  with  B  =  0  and  A  yf  0.  Trajectories  initiated  along  this  line  will  tend  towards  A  — >  oo. 

In  the  next  case,  the  flow  parameters  are  Re  =  5000,  A  =  2.179  x  10 ~5 ,  a  —  100000, 
and  0  =  2.0.  In  Figure  10  pairs  of  trajectories  with  initial  values  of  A  differing  by  0.001  are 
plotted  for  various  initial  values  of  B.  In  many  ways,  the  qualitative  behaviour  seen  here  is 
quite  similar  to  that  illustrated  in  Figure  8.  The  origin  is  an  unstable  node,  there  is  a  saddle 
point  with  nonzero  A  and  B,  and  there  is  a  stable  node  with  nonzero  B  but  A  =  0.  A  semi¬ 
infinite  curve  extends  from  the  origin,  through  the  saddle  point  and  separates  the  domain 
into  a  region  which  is  attracted  to  the  stable  node  with  A  =  0,  B  —  7.7  x  10~4  and  a  region 
in  which  A  —*  oo.  This  shows  how  strong  the  nonlinear  interaction  can  be.  Without  a  Dean 
disturbance,  the  TS  wave  has  unbounded  growth;  however,  the  inclusion  of  an  additional 
disturbance  at  sufficiently  large  amplitude  can  completely  stabilize  the  otherwise  growing 
TS  wave. 

We  conclude  this  section  by  comparing  our  Landau  coefficients  with  those  of  DHZ  for 
two  specific  cases.  Tables  2  and  3  show  the  Landau  constants  from  the  current  theory 
and  from  the  theory  of  DHZ  after  conversion  to  the  current  nondimensionalizations.  The 
quantities  a0io  and  6o,o  from  DHZ  were  taken  directly  from  their  linear  eigenvalue  solver, 
rather  than  by  perturbing  away  from  the  neutral  stability  curve.  Except  for  the  coefficient 
6i (0,  which  represents  the  effect  of  the  TS  wave  on  the  Dean  disturbance,  the  coefficients 
show  satisfactory  agreement.  Singer  and  Zang  [9]  suggested  that  a  yet  unknown  numerical 
problem  in  DHZ  led  to  erroneous  values  of  b  10.  The  large  discrepancy  here  with  respect  to 
bi  o  adds  support  to  that  suggestion. 
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<*0,0 

<*1,0 

<*0,1 

^0,0 

&1,0 

^0,1 

Current 

1.21  x  10"4 

-8.25 

-95000 

3.97  x  10"5 

-482 

-7640 

DHZ 

1.20  x  10“4 

-8.85 

-101000 

3.97  x  10~5 

1.23  x  105 

-8470 

Table  2:  Comparison  of  Landau  Coefficients  for  Re  =  6291.67,  A  =  2.189  x  10  5,  a  —  74257, 
and  (3  —  4.508 


a0,0 

<*1,0 

<*0,1 

&0.0 

&1,0 

^0,1 

Current 

1.18  x  10~5 

3.97 

-46600 

7.66  x  10~6 

-233 

-3340 

DHZ 

1.33  x  10~5 

3.97 

-46900 

7.66  x  10~6 

1.44 

-3390 

Table  3:  Comparison  of  Landau  Coefficients  for  Re  =  4175,  A  =  2.179  x  10  5,  a  —  86800, 
and  j3  =  2.76 

IV.  Comparison  with  Direct  Numerical  Simulation 

We  use  direct  numerical  simulations  (DNS)  of  the  full  Navier-Stokes  equations  to  determine 
the  usefulness  of  the  weakly  nonlinear  theory.  The  numerical  simulation  code  [9]  uses  a  curved 
channel  variant  of  the  method  described  by  Zang  and  Hussaini  [16]  with  the  nonlinear  terms 
in  skew-symmetric  form  [17].  Table  4  reports  the  complex  growth  rates  {a0  0  —  ig0  0)  obtained 
with  a  spectral  linear  stability  code,  and  the  results  of  DNS  with  two  different  time  steps. 
Table  5  provides  details  of  the  parameters  for  the  different  cases.  All  codes  used  65  points 
across  the  channel.  Data  from  the  DNS  were  obtained  after  100  time  steps.  Using  a  time 
step  At  =  0.0001,  the  complex  growth  rates  differ  from  those  predicted  by  linear  theory  by 
less  than  1  part  in  104.  In  the  simulations  reported  below,  we  are  primarily  concerned  with 
long-time  trends  and  final  steady-state  solutions,  hence,  the  additional  time  advancement 
errors  associated  with  the  larger  time  step,  At  =  0.01,  are  not  expected  to  be  important. 
The  larger  time  step  gave  a  maximum  CFL  number  of  0.01,  still  well  within  the  numerically 
stable  regime.  In  one  case  described  below,  the  time  step  was  changed  from  At  =  0.01  to 
At  =  0.0001  after  the  steady  state  solution  was  reached.  There  was  no  evidence  that  the 
solution  was  changing,  even  after  20000  time  steps.  Additional  details  of  code  validation 
studies  in  both  the  linear  and  nonlinear  regimes  are  given  by  Singer  and  Zang  [9]. 


Case 

Linear  Stability 

DNS  At  =  0.01 

DNS  At  -  0.0001 

1 

(3.1189,-4236.9) 

(3.0736,-4241.8) 

(3.1185,-1236.9) 

2 

6.6331 

6.6336 

6.6335 

3 

(  —  138.2126,  —16214.84) 

(-138.2154,-16214.84) 

(-138.2126,-16214.84) 

4 

(86.4287,-2709.11) 

(86.4204,-2709.44) 

(86.4291,-2709.14) 

5 

21.5591 

21.5599 

21.5599 

6 

(-90.61188,  -12993.29) 

(-90.61307,-12993.3) 

(-90.61188,  -12993.34) 

1'able  4:  Comparison  of  Complex  Growth  Rates  (all  values  multiplied  by  104) 
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Case 

Re 

a 

P 

1 

5000 

100000 

2 

5000 

1 

3 

5000 

100000 

2.0 

4 

10000 

80000 

5 

10000 

6 

10000 

80000 

Table  5:  Specification  of  parameters  for  cases  in  Table  4 

Sufficient  spatial  resolution  in  the  simulation  can  be  ensured  using  the  guideline  suggested 
by  Krist  and  Zang  [18].  They  suggest  that  “grid  refinement  is  needed  in  any  direction  when 
the  tail  of  the  energy  spectrum  reaches  10~8  of  the  low-frequency  value.”  This  guarantees 
that  truncation  errors  in  the  velocity  will  be  less  than  0.01%.  Such  detailed  resolution  is 
not  necessary  for  the  purposes  of  these  simulations.  We  found  that  the  tail  of  the  energy 
spectrum  could  be  as  much  as  10-4  of  the  low  frequency  value  and  still  provide  similar  results 
to  those  obtained  when  the  Krist  and  Zang  guideline  was  strictly  followed.  Our  relative 
insensitivity  to  resolution  is  attributed  to  the  rather  minor  role  that  small  scales  play  at  this 
stage  of  transition.  In  what  follows,  using  M  Chebyshev  modes  implies  that  there  are  M  +  1 
Chebyshev  collocation  points  across  the  channel;  using  N  Fourier  modes  implies  that  there 
are  N  collocation  points  in  the  given  direction  so  that  the  Fourier  wavenumbers,  n,  have  a 
range  of  —N/2  +  1  <  nL/( 2tv)  <  N/2  —  1,  where  L  is  the  length  of  the  domain.  In  most 
of  the  simulations  with  Re  =  6291.67,  48  Chebyshev  modes  are  used  to  resolve  the  flow  in 
the  wall-normal  direction,  though  in  the  simulations  with  Re  =  5000,  64  modes  are  used. 
In  the  streamwise  and  spanwise  directions,  2-12  Fourier  modes  are  used,  depending  on  the 
purpose  of  the  given  simulation.  The  code  is  structured  such  that  at  least  2  modes  must 
be  used  in  each  direction.  In  cases  where  no  TS  mode  is  included,  2  modes  are  used  in  the 
streamwise  direction.  Similarly,  where  no  Dean  mode  is  included,  2  modes  are  used  in  the 
spanwise  direction. 

The  initial  disturbances  for  the  direct  numerical  simulations  contain  only  contributions 
from  the  linear  solutions.  The  distortion  of  the  mean  flow  and  higher  harmonic  contributions 
are  not  included  initially;  they  develop  as  the  simulation  progresses. 

Direct  numerical  simulation  is  used  to  study  two  of  the  three  flows  described  in  the 
previous  section.  The  case  in  which  there  is  a  stable  spiral  mode  is  not  simulated  because  the 
time  required  for  a  single  orbit  is  so  large  that  the  computational  cost  would  be  prohibitive. 

We  first  look  at  the  case  with  Re  =  6291.67,  A  =  2.189  x  10-5,  a  =  74257,  and  (3  —  4.51 . 
Four  simulations  were  performed  with  different  initial  conditions;  the  first  contained  only 
a  Dean  vortex,  the  second  contained  only  a  TS  wave,  the  last  two  contained  both  types  of 
disturbances. 

When  only  a  Dean  disturbance  is  included  in  the  initial  conditions,  the  weakly  nonlinear 
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theory  suggests  that  an  equilibrium  state  with  B  =  7.0  x  10-5  develops.  This  corresponds 
to  a  maximum  streamwise  velocity  perturbation  of  5.84%  of  U*.  We  initiated  the  simulation 
with  the  Dean  eigenfunction  obtained  from  linear  stability  theory  with  an  initial  strength, 
B  =  6.95  x  10-5.  Twelve  Fourier  modes  were  used  in  the  spanwise  direction,  though  half 
the  modes  had  energies  which  were  more  than  8  orders  of  magnitude  less  than  the  primary 
mode.  These  higher  order  modes  had  a  negligible  effect  on  the  simulation.  At  the  start  of 
the  simulation,  the  instantaneous  growth  rate  of  the  Dean  disturbance  was  approximately 
equal  to  its  linear  value.  At  the  time,  t  =  34845,  the  magnitude  of  the  instantaneous  growth 
rate  had  decreased  by  more  than  a  factor  of  100  and  the  simulation  was  stopped  with 
B  —  6.7  x  10~5.  In  this  case,  the  difference  between  the  equilibrium  state  predicted  by  the 
weakly  nonlinear  theory  and  that  obtained  from  the  DNS  is  less  that  5%  and  is  considered 
quite  good. 

With  only  a  TS  wave,  we  expect  to  obtain  a  TS  equilibrium  state  with  A  —  3.8  X  10-3, 
resulting  in  a  maximum  streamwise  disturbance  of  1.94%  of  U*.  The  simulation  was  initiated 
with  A  =  3.5  x  10-3  and  was  stopped  at  a  time  t  =  11367.  At  this  time,  A  =  4.0  x  10-3. 
Figure  11  illustrates  the  time  evolution  of  the  amplitude  of  the  wave.  The  dashed  line 
shows  the  equilibrium  amplitude  predicted  by  our  weakly  nonlinear  theory.  Note  that  as  the 
amplitude  approaches  its  equilibrium  value,  the  growth  rate  decreases.  Over  480  periods  of 
the  TS  wave  were  simulated  for  the  data  presented  in  Figure  11.  The  difference  between 
the  equilibrium  amplitude  predicted  by  the  weakly  nonlinear  theory  and  that  obtain  by  the 
DNS  is  about  10%,  which  we  considered  acceptable  agreement.  The  12  Fourier  modes  used 
in  the  streamwise  direction  were  sufficient  for  an  8  order  of  magnitude  drop  in  the  energy. 

We  performed  two  simulations  with  both  Dean  and  TS  perturbations  included  in  the 
initial  conditions.  In  these,  we  used  6  modes  in  both  the  streamwise  and  spanwise  directions 
and  48  Chebyshev  modes  in  the  wall-normal  direction.  In  these  simulations,  the  spanwise 
resolution  was  only  sufficient  to  allow  an  energy  drop  of  105,  wdiile  the  streamwise  resolution 
allowed  as  little  as  a  104  drop  in  the  energy.  This  would  lead  us  to  suspect  errors  in  the 
velocity  on  the  order  of  1%.  To  determine  whether  this  could  lead  to  erroneous  physical 
results,  we  re-simulated  our  case  which  had  the  largest  velocity  perturbations  using  12  modes 
in  the  streamwise  and  spanwise  directions  and  64  modes  in  the  wall-normal  direction.  The 
er  ergy  spectrum  in  the  highly  resolved  case  exhibited  an  energy  drop  of  1 09  in  the  streamwise 
direction  and  1014  in  the  spanwise  direction.  The  evolution  of  the  disturbance  amplitudes 
was  the  same  as  that  observed  in  the  less  resolved  case.  In  both  cases  we  initiated  the  Dean 
disturbances  with  B  =  6.95  x  10-5.  This  value  is  slightly  less  than  its  equilibrium  amplitude 
which  we  obtained  above. 

In  the  first  case,  the  initial  strength  of  the  TS  disturbance  was  A  —  2.94  x  10-4,  approx¬ 
imately  1/13  its  predicted  equilibrium  value.  The  long  time  behavior  of  this  case  required 
12  Fourier  modes  in  the  spanwise  direction.  The  amplitude  evolution  of  the  disturbances  is 
illustrated  in  Figure  12a.  Here  we  plot  the  amplitude  of  both  the  Dean  and  TS  disturbances, 
normalized  by  their  respective  amplitudes  at  time  t  —  0.  The  Dean  disturbance  develops  to¬ 
wards  its  equilibrium  state  while  the  TS  disturbance  decays.  In  Figure  12b  the  predictions  of 
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the  weakly  nonlinear  theory  are  plotted.  The  agreement  between  the  theoretical  predictions 
and  the  DNS  results  is  very  good. 

In  the  second  case,  the  initial  strength  of  the  TS  disturbance  was  A  =  2.94  x  1 0 “ 3 ,  ap¬ 
proximately  10/13  its  predicted  equilibrium  value.  Figure  13a  shows  the  amplitude  histories 
for  this  case.  Unlike  the  previous  situation,  here  it  is  the  TS  wave  which  grows  towards  its 
equilibrium  value  and  the  Dean  disturbance  which  experiences  rapid  decay.  The  predictions 
of  our  weakly  nonlinear  theory  are  plotted  in  Figure  13b.  Again  the  agreement  between 
theory  and  simulation  is  quite  good. 

The  next  flow  which  we  investigated  with  numerical  simulations  corresponds  to  the  last 
case  studied  in  the  previous  section.  Here  Re  =  5000,  A  =  Ac  =  2.179  x  10~5,  a  =  100000 
and  (3  =  2.0.  Simulations  were  performed  with  a  TS  wave  only,  a  Dean  vortex  only,  and  a 
combination  of  both  disturbances. 

The  weakly  nonlinear  theory  predicts  that  a  TS  wave  alone  in  this  flow  becomes  more 
unstable  as  its  amplitude  increases  since  both  ao,o  and  a^o  are  positive.  A  numerical  simu¬ 
lation  initiated  with  A  =  0.002  and  B  =  0  confirmed  that  the  amplitude  A  increased  with 
increasing  growth  rate,  at  least  until  A  =  0.0048  when  the  simulation  was  terminated.  Six 
Fourier  modes  were  used  in  the  streamwise  direction. 

A  Dean  vortex  alone  in  this  flow  is  predicted  to  reach  an  equilibrium  state  with  B  =  7.7  x 
10-4.  This  value  of  B  corresponds  to  a  maximum  streamwise  velocity  perturbation  which  is 
62.4%  of  the  undisturbed  laminar  bulk  velocity.  An  equilibrium  state  with  B  =  4.3  x  10~4 
was  obtained  in  the  direct  numerical  simulations.  Twelve  Fourier  modes  were  used  in  the 
spanwise  direction.  The  difference  between  the  DNS  result  and  that  of  the  weakly  nonlinear 
theory  is  not  unexpected  since  the  disturbance  amplitudes  are  so  large. 

One  final  simulation  was  performed  with  the  initial  values  of  A  and  B  chosen  as  A  =  0.002 
and  B  =  0.00065.  In  this  case,  the  weakly  nonlinear  theory  predicts  that  the  TS  wave  will 
decay,  leaving  only  a  Dean  disturbance  in  the  flow.  The  time  evolutions  of  the  TS  and  Dean 
disturbances  are  plotted  in  Figure  14a  through  t  —  500.  Six  Fourier  modes  were  used  in 
the  streamwise  direction,  twelve  in  the  spanwise.  At  this  time  the  spectral  energy  decay 
in  the  spanwise  modes  is  105.  Though  the  simulation  was  terminated  before  a  steady-state 
solution  was  reached,  the  prediction  of  the  weakly  nonlinear  theory  (plotted  in  Figure  14b  is 
qualitatively  supported.  After  the  initial  transient,  the  TS  wave  decays  rapidly.  Considering 
the  large  amplitude  disturbances  in  the  flow,  the  theory  and  the  simulation  agree  well. 


V.  Conclusions 


We  have  presented  a  weakly  nonlinear  theory  to  describe  interactions  between  TS  waves  and 
Dean  vortices  in  curved  channel  flow.  The  approach  used  to  calculate  the  Landau  coefficients 
is  an  extension  of  that  developed  by  Herbert  [10,  11].  This  approach  provides  not  only 
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the  Landau  coefficients,  but  also  the  shape  change  of  the  fundamental  disturbances  with 
amplitude.  A  standard  Chebyshev  collocation  procedure  is  used  to  obtain  numerical  results. 
We  have  restricted  our  study  to  very  weakly  curved  channels,  where  the  critical  Reynolds 
numbers  for  both  TS  and  Dean  disturbances  are  approximately  equal.  For  low  Reynolds 
numbers  the  interaction  effects  are  generally  stabilizing  to  both  types  of  disturbances.  In 
the  case  of  high  Reynolds  numbers  the  theory  predicts  that  the  interaction  can  destabilize 
the  TS  wave.  Destabilization  of  the  Dean  vortex  due  to  the  TS  wave  occurs  only  for  very 
low  wavenumber  Dean  disturbances. 

Three  different  scenarios  of  TS-Dean  interactions  have  been  studied  in  detail.  Phase 
plane  trajectories  indicate  the  existence  of  stable  and  unstable  critical  points.  Two  cases 
exhibit  saddle  point  behaviour  while  the  third  has  a  stable  spiral  node. 

Direct  numerical  simulations  verify  the  predictions  of  the  theory.  When  the  amplitudes  of 
the  disturbances  are  not  too  large,  the  direct  numerical  simulations  and  the  weakly  nonlinear 
theory  agree  quite  well.  When  the  amplitudes  of  the  disturbances  are  large  (on  the  order  of 
60%  of  the  centerline  velocity),  the  simulations  and  the  theory  still  qualitatively  agree. 

Comparison  of  our  results  with  those  of  DHZ,  who  also  developed  a  similar  weakly  non¬ 
linear  theory  strongly  suggest  that  DHZ  have  a  numerical  error  in  the  computation  of  one 
of  their  Landau  coefficients. 
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Appendix  A 


APPENDICES 


I.  The  Left  Hand  Side  Operator  Ln/ 


Ln,Z  — 


/L1,1  L1'2  Z1,3  Z1,4  \ 

L2'1  L2'2  L2'2  Z2'4 
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II.  The  Right  Hand  Side  Vectors  Rs 


R,= 


Equation  system  5  corresponds  to  the  TS  harmonic.  Its  right  hand  side  vector  is  given 


R\  =  —iaug1  ’0>0i0ur  1  ’°'0'0/ r 

+  (ufl1'0’0'0)2/r 

d  ,  ~  i,o,o,o\  -  1,0, 0,0 

-—{Ur  )ur 

dr 

R\  =  -za(uV'0'°’0)2/r 

-uel'0'0'0Urlfl'0fllr 
_A(uy.°.  O.OJ^I.O.O.O 

Rl  =  —  ZQUfl1'0,0’0uz1,0,0'0/r 


+  n'r1, 0.0.0  “  (tf  l’0,0'0) 
dr 


Ri  =  0 


Equation  system  6  corresponds  to  the  Dean  harmonic.  Its  right  hand  side  vector  is  given 


-ipuz0fi'lfiur0'0X0 
+  (^°'°'1'°)2/r 

d  ,  *  o,o,i,o\  »  o,o, i,o 

~Tr(Ur  K 

•a  -  0,0,1, 0,«  0,0, 1,0 
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dr 
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.  nnin  d  ,  ,  noin. 


0,0, 1,0 


f(Wz0,0,1,°) 


Equation  system  7  corresponds  to  the  product  of  the  TS  and  Dean  fundamentals.  Its 
right  hand  side  vector  is  given  by: 

r>  1  •  A  0,0. 1,0  *  1, 0,0,0  / 

Ii7  —  —  iqlxlq  ’  '  '  ur  *  ‘  ’  /r 


22 
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_{d  o  +  d  ^ 

dr  dr 

— ia  tig  0,0,1 ,0tlz  1  •0,0,0/r 

-i/Ju,0*0-1-0^1-0*0*0 

_(i-(^o.o.i,o)u;llo.o.o  +  A(^1.0.0,0)u%0,0,1.0) 
dr  dr 


Equation  system  8  corresponds  to  the  product  of  the  TS  fundamental  and  the  complex 
conjugate  of  the  Dean  fundamental.  Its  right  hand  side  vector  is  given  by: 


—ia  rig0,0'-1 ,0ur  1 ,0,0,0/r 
+f/3ur0,0’_1,0uI1, 0,0,0 
+2  de0,0,-1,%1,0,0,0/r 

ar  ar 


=  —iaug 


O.O.-l.O l.O.O.O 


+  2/Sufl0,0,-1,0UI1, 0,0,0 

—  (ug0,0’-1,0^1, 0,0,0  +  u91,0,0,0ur0,0,-1,0)/r 


=  —iaug 


0,0, -1,0"  1 ,0,0,0. 


+  i/?'!iz0,0,-1,0Uz1, 0,0,0 
_(iL(t2,o.o.-i.o)tiri.o.o.o  +  l(u*^.o.o.o)uVo,o,-i,o) 


Equation  system  9  corresponds  to  the  self  correction  of  Ere  TS  fundamental.  Its  right 
hand  side  vector  is  given  by: 


-za(2u9-1,0,0,0ur2, 0,0,0  +  ugu,1,u,V1,u,u,u  -  Uy,u‘ 

+  (2us-1,0,0,%2, 0,0,0  +  2  ug0,1,0,0u91,0,0,0)/r 

_(d(u-r".».o.o)tfr-i.o.o.o  + 

dr  dr 


0, 1,0,0  -  1,0, 0,0  -  2, 0,0,0  7  -1,0,0  0W 

9  lir  U>Q  dr  J  /  7 
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Rl  = 


Rl 


R$ 


+  ±(ur0‘1’0fi)uT1'0fi’°  +  y-(t2r1,0,0’0)Ur0,1’0,0) 
dr  dr 

-(ai>0  -  igi.o)^1’0,0'0 

—ia(ug~l’0'0’0ug2,0'0'0  +  ug0’l-0'W'0’0’0)/r 

-(ue-1A0’°Ur2‘0A0  +  ue3'0'0'0ur~1>0’0’0  +  Ug0'1'0'0^1’0’0'0  +  US— Ur 


0,1, 0,0,;  1,0, 0,0  J_  ,7,1.0, 0.0,-  0,1,0, Oyr 


2,0,0,0s  «  -1,0, 0,0 


o)tira. 0,0,0  +  ^.(tI/,o,o,o)t2r 

N  dv  ri'r 


dr 


+  — -(t2fl0,1,0,0)ttI.i’0'0,o  +  ~  (d*1'0'0'0)^0'1'0'0) 
dr  dr 

/  \  -  1,0, 0,0 

—  (^1,0  —  151,0  jag 

=  — za(2ue_1,0’0,0tt*2’0,0,0  +  -  t4wClM°)/r 

-(i(«*"1,0,0,0K2*' O'0,0  +  ^(uz2fi-ofi)ur-w 


dr 


+  —(uz°’l'0fi)urlfi‘0fi  +  f-(uVl0,0-0)Ur0ll’0l°) 

dr  dr 

—  (ai,o  —  151,0)^ 

0 


1,0, 0,0 


Equation  system  10  corresponds  to  the  self  correction  of  the  Dean  fundamental.  Its  right 
hand  side  vector  is  given  by: 


*10 


=  —i/3(2uz0'0'~1'0ur0'0'2'0  +  uz°’0'0'luT°>0’1’0  -  u^’W’0’-1’0) 
+  (2ue0'0’-1'V'0'2’0  +  2ue°-O'O-1ii0O'O'1'°)/r 


R2 


10 


R3 


10 


Ue 

-( Afu,0'0’2'0)*0'0’-1'0  +  ^-(«r0'°'-'-0)ur0'0'2'0 

ar  dr 

+  ±(ur0-°’°'l)ur  O'O’I’O  +  ^(u,0-0-1-0)^0,0,1) 

dr  ar 

/,  •  >  \  -  0,0, 1,0 

—  (.00,1  —  lao.lJUr 

-z/?(2u90'0'2'0^0'0'-1'0  +  U*0A1' V’0’0'1  -  tf,0-0-'1-0^0-0-2-0) 

-(dfl°'°'-1-V0'0'2'0  +  ug0,0’2’0ur0,0’-1’0 

+U90,0’0,1ur0,0'1'0  +  ug0,0’1’0ur0’0,0’1)/r 

f±^B°'0'-l'°)ur0'  °’2*0  +  A(tf,°'0'2'0)ur°'0'-3'0 

' dr '  dr 

+  — (uV- °l01 1)0  °-1>0  +  ^-{ug0'0’1'0)^0'0'0'1) 

-(t,  -  1*0,1  Ms0'0'1’0 

-^(ti,0-0*-1*0!*,0-0-210  +  «,0'0,0,l«.0’0,1’0) 

-(  A(di0’0--1’0)dr0'0'2'0  +  ^-(u,0-0-2-0)ur0-0*-ll° 

v  dr  dr 

,  d  (  -  0,0,0,1\  -  0,0, 1,0  ,  ^  /  7  0,0,1,0s  7  0,0,0,1s 

+*:(“'  )u'  +*(u-  ,Ur  J 


/ »  •  i  \  A  0,0,1 ,0 

~(oo,l  —  1*0,1  )uz 
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R,  o  =  o 


Equation  system  11  corresponds  to  the  correction  of  the  TS  fundamental  by  the  Dean. 
Its  right  hand  side  vector  is  given  by: 

Rh  =  — ta(u80’0,_1,0txP1,0,1*0  +  us0,0,0,1^1-0,0,0  +  titf0.0.i.0t2pi.o.-i.0)/r 
4-i/3(ur0,0,-1,0uz1, 0,1,0  _  wV0'0'1'0^1'0'-!'0 
+dr1,0,_1,(K0, 0,1,0  -  ^1,0.1, 0^0,0,- 1.0) 

+2(u80,0’-1,0Ufl1, 0,1,0  +  u90,0,0,1ua1,0,0'0  4-  ue°'0,1,0t2e1,0'-1,0)/r 
-(^(u;0,0'-1,0)^1,0,1,0  +  i.(u;o,o,o,i)tiri, o.o.o 

+  ^-(ur0,0,1,0)ur1,0,-1,°  +  i(tir1'0,-1'°)ur0, 0,1,0 

+  i_(^l,0,0,0)^0, 0,0,1  +  iL(^1.0.1. 0^.  0,0, -1,0s 

dr  dr  ' 

“(«0,1  -  *00,1  )ur1,0'0,0 


0,0,_1,0Ufl1, 0,1,0  +  ti*0,0,0,1^1, 0,0,0  -t-  i;fl0,0,1,0tifl1'°-i'0)/r 


P2  -  - llu|l|u  I  *  u,U,U,i  ~  l.U.U.U  , 

riu  —  layue  Ug  4-  Ug  '  Ug  '  '  '  +  Ug 


4-tjS(ti90,0--,,V,( 5,1,0  -  u80,0,1,0uz1,0,-1,° 
W’°’~1,0u20, 0,1,0  -  tzfl1,0,1,0u20,0,~1,0) 


•  ,0,0  0 


4-ite0,0,1,0ur1,0,-1,°  f  V'°,_1'V0, 0,1,0 


4 -uV,0,0,0u;0, 0,0,1  4-  «81,0,1,0ur0,0,_1,0)/i 


d_ 

'dr 

A 

dr 

d 


dr 

d 


-(u80,0,0,1)ur1,0,0'° 


+  —  (ua0,0,1,0)^1,0'-1,0  4-  ^(tx81,0,“1,0)ur0, 0,1,0 
»?„i>°>°,0^;  o.o.o, i  _l  Af.:  i.o.i, 0\„  ?  o,o,-i,o^ 


4--(t281,0,0,0)ur0,0,0,i  +  f  (u81,0,1,0)ur0,0,-1,°) 

-(«o.i  -  *<7o,iK1,0,0,0 

Rh  =  -ia(u80,0,-1,0t2z1,0,1,0  4-  u*0,0,0,1^1, 0,0,0  +  ue°.°.1,0u21,0,-1,0)/r 

-(^(^0'0,“1,0)Ur1, 0,1,0  4  ^(U,0,0,0,1)^1, 0,0,0 

+  iu%o,,,)t2r,o-,o  +  ^.(lii1.0.-1.0jlIr0.0.1,0 


dr 

d 


+  -(.V,0,0,0)ur0, 0,0,1  +  ^(txV,0,1,0K0,0.-1.°) 


-(<*0.1  -  470,1k1, 0,0,0 


A4 


li 


0 


Equation  system  12  corresponds  to  the  correction  of  the  Dean  fundamental  by  the  TS. 
Its  right  hand  side  vector  is  given  by: 


Rn  = 


-ia(ur1,0,0,V1,( 3,1,0  +  uV1,0,1,0^1, 0,0,0 


ue 


1 ,0,0,0  ~  -1,0, 1,0  ?^1, 0.1,0,;  -1,0, 0,0 


4-  Ug 


UT 
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r\2 


R 3 


7  -1,0, 1,0,-  1,0, 0,0  J_  ,7  0,0, 1,0^  0,1, 0,0  _L  ,7  1,0, 1,0,7  -1.0, 0,0 


+  Ur 


uz 


-1,0, 1,0,7  1,0, 0,0  ,  ,7  0,0, 1,0,7  0,1, 0,0 


+  Ug 


ug 


-t/3(  — .«  +  ur 

+2('Ufl_1,0'0,0Ufl1, 0,1,0  +  we  - us 

_(i_(Tir-^o.o.o)liri.o,i.o  +  i_(u;-i.o.i,o)u%i.o,o,o 

dr x  ar 

+  -j-(  ur  0,0,1 ,0)ur0,i  ,0,°  +  ^-(rir0,1,0,°)i2r0, 0,1,0 
dr  dr 

+  1.(72/, °.0, 0)^-1, 0,1,0  +  A(t2r1,0,1*0)i2r"1,0,0,°) 

ar  dr 

-(6i,o  -  ^i,o)ur0, 0,1,0 

-i/(uV1|0ll’VA0'°  +  ug0, 0,1,0  u20, 1,0,0  +  ufi1,0,1,0ix2-1,0,0,0) 


-(w« 


7  -1,0, 0,0,7  1,0, 1,0 


Ur 


+  Ug 


i,o, 1,0,7  1,0, 0,0 
ur 


1,7  0,0, 1,0,7  0,1, 0,0  ,  ,7  0,1, 0,0,7  0,0, 1,0 
+  Ug  Ur  +  Ug  ’  Ur 

+«y“°ur-iA1-0 + ri#>.».i.ot!r-..°.o.O)/r 

dr  dr 

+  ±(usw'0)ur°'1M  +  d{us°-'fifi)urw-° 
dr  dr 

+4-(us'A0'°)ur-'M'° + -/uyA,'>K','w) 

dr  dr 

-(Vo  -  i/>1,o)u»'W'° 

1,0, 0,0,7  1,0, 1,0 


12  —  lO(\Ug 

~lP(Uz 


U, 


+  Ug 


-1, 0,1,0^  1,0, 0,0  _  ,7+ ,0,0,0, 7  -1, 0,1,0  _  ,7,1, 0,1, 0,7  -1,0, 0,0 


-  Ug 


U 


-  Ug 


Ur_ 


°)/r 


1,0, 1,0^  1,0, 0,0  _j_  ,7  0,0,1 ,0,7  0,1, 0,0  ,7  1,0, 1,0,7  -1,0, 0,0 


+  u 


u 


+  Uz 


uz 


7?4 

/l12 


ar  ar 

+  i.(u-_l.o,o.o^-l,o,,,o  +  4  (u%  Wu-IMO) 
ar  ar 

-(*1.0  -  *>H,oKOA1’0 
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Figure  6 
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Figure  6:  Zero  contour  of  o0li,  the  effect  of  the  Dean  disturbance  on  theTS  disturbance.  The 
rectangles  represent  the  range  of  a  and  (3  over  which  individual  TS  and  Dean  disturbances 
are  linearly  unstable  at  the  given  Reynolds  number. 
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Figure  7 


a 

Figure  7:  Zero  contour  of  &ii0,  the  effect  of  the  TS  disturbance  on  the  Dean  disturbance.  The 
rectangles  represent  the  range  of  a  and  0  over  which  individual  TS  and  Dean  disturbances 
are  linearly  unstable  at  the  given  Reynolds  number. 


Figure  10:  Trajectories  in  phase-space  with  Re  =  5000,  A  =  2.179  x  10~ 5 ,  a  =  100000,  and 
(3  =  2.0.  Squares  represent  initial  conditions,  filled  circles  represent  stable  critical  points. 
Open  circles  represent  other  critical  points. 


Figure  12(a) 


Figure  12(b) 


time 

Figure  12:  Time  evolution  of  TS( - )  and  Dean( - )  disturbances  for  Re  =  6291.67, 

A  =  2.189  x  10-5,  a  =  74257,  and  0  =  4.51.  The  initial  TS  disturbance  has  A  =  2.7  x  10-,V 
The  initial  Dean  disturbance  has  B  =  7.0  x  10-5.  The  results  from  the  direct  numerical 
simulation  are  shown  in  (a).  The  predictions  of  the  weakly  nonlinear  theory  are  given  in  (b). 
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Figure  13(a) 
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Figure  13(b) 
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Figure  13:  Time  evolution  of  TS( - )  and  Dean( - )  disturbances  for  Re  =  6291.67, 

A  =  2.189  x  10"*,  a  —  74257,  and  0  =  4.51.  The  initial  TS  disturbance  has  A  =  2.7  x  10-3. 
The  initial  Dean  disturbance  has  3  =  7.0  x  10-5.  The  results  from  the  direct  numerical 
simulation  are  shown  in  (a).  The  predictions  of  the  weakly  nonlinear  theory  are  given  in  (b). 
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Figure  14:  Time  evolution  of  TS( - )  and  Dean( -  )  disturbances  for  Re  =  5000, 

X  =  2.179  x  10-5,  a  =  100000,  and  /3  =  2.0.  The  initial  TS  disturbance  has  A  =  2.0  x  10  3. 
The  initial  Dean  disturbance  has  B  =  6.5  x  10“4.  The  results  from  the  direct  numerical 
simulation  are  shown  in  (a).  The  predictions  of  the  weakly  nonlinear  theory  are  given  in  (b). 
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